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Mapping out the quark structure of hadrons in QCD* 
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In the context of QCD sum rules with non-local condensates we present a pion dis- 
tribution amplitude, which is double-humped with its end-points x — > (0 , 1) strongly 
suppressed, and show that it matches the CLEO experimental data on the pion- 
photon transition at the la level accuracy, being also in compliance with the CELLO 
data. We also include some comments on the nucleon distribution amplitude and 
q ' the nucleon evolution equation. 
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I. A TRIBUTE TO PROF. EFREMOV'S CELEBRATION OF HIS 70th 

BIRTHDAY 

Prof. Efremov gives us the opportunity to point out in this Festschrift the influence of his 
work on our own research activities. 

A. V. Efremov is one of the inventors of factorization theorems in quantum field theory 
^ . that are particularly indispensable in applying perturbative QCD in inclusive 0] and exclusive 
reactions 0, HI involving hadrons. Without these tools, the experimental verification of QCD 
would constitute an intractable task. Together with his then student A. V. Radyushkin he 
accomplished the factorization theorems for the meson form factors, linking diagrammatic 
techniques with the operator product expansion (OPE). The grounds for these works were 
supplied by previous investigations by Efremov 0| and Efremov and collaborators 

Moreover, Efremov and Radyushkin have diagonalized the anomalous dimensions matrix 
■ for meson operators (in leading order) in terms of Gegenbauer polynomials and first obtained 
the asymptotic distribution amplitude (DA) ip(x,fi 2 — > 00) — > ip as \x) = 6x(l — x) 0,0- 

Factorization theorems make it possible to calculate various hard processes in QCD 

involving mesons, in which the meson DAs enter as the central nonperturbative input. 

In the context of the present occasion, we are primarily interested in presenting recent 
achievements in describing the pion characteristics by mapping out its internal quark structure. 
A short note on the nucleon is also included. 



II. NON-LOCAL CONDENSATES AND PION DISTRIBUTION AMPLITUDE 

The pion DA of twist-2, <p n (x, /i 2 ), is a gauge- and process-independent characteristic of the 
pion that universally specifies the longitudinal momentum xP distribution of valence quarks 
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in the pion with momentum P 

(0 | d(0)7^ 6 S(0, z)u{z) | tt(P)) = z/ w P" f dxe lx(zP ^ <p v (x, /i 2 ) (2.1) 

and where £7(0, z) = V exp^—igj^ t a A'f l (y)dy^ is a phase factor, path-ordered along the 
straight line connecting the points and z to preserve gauge invariance. 

A. Average QCD vacuum quark virtuality Aq 

The pion DA encapsulates the long-distance effects and therefore reflects the nonpertur- 
bative features of the QCD vacuum. The latter can be effectively parameterized in terms of 
non-local condensates, as developed in [H, 0, by A. Radyushkin and two of us (A.B. and 
S.M.). This provides a reliable method to construct hadron DAs that inherently accounts for 
the fact that quarks and gluons can flow through the QCD vacuum with non-zero momen- 



tum k q . This means, in particular, that the average virtuality of vacuum quarks, (A; 2 ) = A 2 



is not zero, like in the local sum-rule approach |llj |. but can have values in the range [1 
A2 = (q(iga tiU G tlu )q)/(2(qq)) = 0.35 - 0.55 GeV 2 . Therefore, the non-local condensates in 
the coordinate representation, say, (q(0)E(0, z)q(z)), are no longer constants, but depend on 
the interval z 2 in Euclidean space and decay with the correlation length A ~ 1/A g . Lacking an 
exact knowledge of non-local condensates of higher dimensionality, one has de facto to resort 
to specific Anzatze ^3|, in order to parameterize the non-local condensates. Nevertheless, it 
is important to stress that we were able to determine in 0| A 2 directly from the CLEO data 
lif within the range predicted by QCD sum rules 12|| and lattice simulations [13], favoring 



the value A 2 ~ 0.4 GeV 2 . 

B. QCD sum rules 

The distribution amplitudes ^(^(x, /i 2 ) for the pion and its first resonance can be related 
to the non-local condensates by means of the following sum rule that is based on the correlator 
of two axial currents 

ffrM + /W*) exp{-^} = [* P pert (x; s)e-^ l2 ds + $ G (x; M 2 ) 

J 

+ ai jur? E ^ ) . ( 2 - 2 ) 



i=S,V,T 1 ,T 2 ,T 3 



where the index i runs over scalar, vector, and tensor condensates 0, 17], M 2 is the Borel 
parameter, and s° is the duality interval in the axial channel. Above, the dependence on the 
non-locality parameter enters on the RHS in the way exemplified by the numerically important 
scalar-condensate contribution 

18 ( 

$ 5 (x; M 2 ) = -£-^2\9 (.x > A > x) x [x + (A - x) In (x)) + (x x) + 

+6(1 > A)6 (A > x > A) [A + (A - 2xx) ln(A)] } (2.3) 

with A = A 2 /(2M 2 ), A = 1 — A and x = 1 — x. In the so-called local approach [11], the 
end-point contributions (x —> or 1) are strongly enhanced by S(x),5'(x) . . . because they 
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disregard the finiteness of the vacuum correlation length A by setting in Eq. ()2.3|) A 2 , — > to 
obtain 

hm $ 5 (x; M 2 ) = 9 [5{x) + 5(1 - x)} . (2.4) 

In contrast, taking into account the non- locality of the condensates via A 2 ,, leads to a strong 
suppression of these regions. Due to the end-point suppression property, the sum rule ()2.2|) 
allows us to determine the first ten moments (£ N } W = Jq (p 7r (x)(2x — 1) N dx of the pion DA 

and independently also the inverse moment (x^ 1 )^ = J Q <p n (x) x~~ x dx quite accurately (see in 
Til igj | for more details). The intrinsic accuracy of this procedure admits to obtain the pion 



DA moments with uncertainties varying in the range of 10%. 



C. Models for the pion distribution amplitude 
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Models for the pion DA, in correspondence to the extracted moments, can be constructed 

in different ways [1, [13]. However, on 
the grounds explained above, it appears 
that two-parameter models, the param- 
eters being the first Gegenbauer coeffi- 
cients a 2 and d4, enable one to fit all the 
moment constraints for (£ N )n, as well as 
to reproduce the value of {x' 1 )^ within 
the QCD sum-rule error range, resulting 
into a "bunch" of DAs displayed in Fig. ^ 
The optimum sample out of this "bunch" , 
termed BMS model, is described by the 
following expression 




FIG. 1: Comparison of selected pion DAs denoted by ob- 
vious acronyms: tp as (dotted line) 0.0- Vpr (dashed line) 
[iof , t^Dor (dot-dashed line) , and </?bms (solid line) Q . 
Also shown is the whole "bunch" determined via QCD sum 
rules with non-local condensates Q|. All DAs are normal- 
ized at the same scale /l*q » 1 GeV 2 . 
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+ a 4 cf 2 (2z 



l + a 2 G 2 3/2 (2x- 1) 
I i (2.5) 



with a 2 = +0.20, a 4 = —0.14 and is em- 
phasized by a solid line in Fig. ^ The 
shape of this "bunch" is confirmed by a 
non-diagonal correlator, based on the QCD sum rules considered in 2^ . 



III. CLEO DATA ANALYSIS 



The CLEO data 2| on provide one rigorous constraint on theoretical models for the 



7T7 

pion DA in QCD. Indeed, it was first shown in |23j that these data exclude the CZ pion DA 
because the prediction derived from it overshoots these data by orders of magnitude. Very 
recently, we analyzed [lil Il9l | the CLEO data by combining attributes from QCD light-cone 
sum rules jiil |25| . NLO Efremov-Radyushkin-Brodsky-Lepage (ERBL) 0, H E| evolution 
26l 27], and detailed estimates of uncertainties owing to higher- twist contributions and NNLO 
perturbative corrections [28j]. 

The upshot of this analysis is that the CZ pion DA is excluded at the 4a level 
of accuracy and — perhaps somewhat surprisingly — that also the asymptotic pion DA 
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lies outside the 3a error ellipse in the 
(ci2, 04) plane (see Fig. |2J), even if one al- 
lows the theoretical uncertainties owing to 
unknown higher-twist contributions to be 
of the order of 30% and presumes that the 
size of NNLO perturbative corrections is 
also large. On the other hand, the BMS 
pion DA calculated with a vacuum virtu- 
ally A 2 , ~ 0.4 GeV 2 was found to be in- 
side the la error ellipse, while other rival 
models, based on differing instantons ap- 
proaches 0, EH , or derived with the aid 
of lattice simulations [3(|, are located in 
the vicinity of the border of the 3a con- 
tour. It is worth emphasizing that the 
more precise the instanton-based models 
become, the further away from the asymp- 
totic pion DA towards the region of the 
"bunch" they move (we refer to \lj^ for 
more details). 1 It was pointed out be- 
fore in [32] that the CLEO data ask for 
a broader pion DA than the asymptotic 
one. 

In Fig. El (left panel) we compare 
our prediction for the scaled pion-photon 
transition form factor with those from the CZ model (upper dashed line) and the asymptotic 
DA (lower dashed line). One observes that the strip obtained from the "bunch" of DAs is 
in very good agreement with both the CLEO data and also with the CELLO data [33]. The 
right panel of this figure illustrates in the form of a shaded band the region of uncertainty 
induced by our limited knowledge of higher-twist contributions. One observes that even the 
low-Q 2 CELLO data are in reasonable compliance with the theoretical prediction (the shaded 
strip). 

Let us close this section by mentioning that other _approaches claim to be able to describe 
the CLEO data with the asymptotic pion DA (23I. I341 135^ . taking into account only the leading- 
twist contribution and using only perturbative QCD (see for more details in |14j|). 

IV. OTHER EXCLUSIVE PROCESSES 

Factorization theorems can be extended — at least formally — to baryons and their form 
factors ja|. The primary goal below is to give a brief summary of main results rather than 
to review the subject and the status of individual exclusive processes or baryon DAs (for a 
recent review, we refer to j2f3|). For instance, the situation concerning the nucleon DA is more 



The new model relative to |3l| . proposed in |2l| . involves more than two Gegenbauer coefficients and can 
therefore not be displayed in Fig. [2] However, reverting this model to an approximate one by utilizing only 
two (effective) Gegenbauer coefficients 02 and 04 shows that it is close to the 3a error ellipse boundary, as 
said above. 
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FIG. 2: Analysis of the CLEO data on F 7!1 * 1 {Q 2 ) in the 
(0,2, 04) plane in terms of error regions around the best-fit 
point (blue cross) with the following designations: la (bro- 
ken green line); 2a (solid blue line); 3a (dashed-dotted red 
line). Various theoretical models are also shown for com- 
parison. The designations are as follows: ♦-the asymptotic 
DAjX-BMS model, B-CZ DA, +-best-fit point, it 0, 
♦ [20L A j^-instanton models, and T-transverse lattice 
result |30| . The slanted green rectangle represents the BMS 
"bunch" of pion DAs dictated by the nonlocal QCD sum 
rules for the value A 2 = 0.4 GeV 2 . All constraints are eval- 
uated at fi 2 = 5.76 GeV 2 after NLO ERBL evolution. 
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FIG. 3: Left: Light-cone sum-r ule p redictions for Q 2 F 7 * 1 ^ 7T (Q 2 ) in comparison with the CELLO (diamonds, 
|3^ 1 and the CLEO (triangles, Q|) experimental data evaluated with the twist-4 parameter value S 2 „ A = 
0.19 GeV 2 ElEi- 

The predictions correspond to selected pion DAs; viz., (fez (upper dashed line) 
BMS- "bunch" (shaded strip) Q, and tpas (lower dashed line) 0,0]. Right: Our prediction for Q 2 F r >~ f _ +7r {Q 2 ) 
corresponding to the "bunch" of pion DAs in Fig. ^ (shaded strip) in comparison with experimental data for 
twist-4 parameter values varied in the range i£f w _ 4 = 0-15 — 0.23 GeV 2 . 



controversial compared to the meson case. It is undoubtedly true that the asymptotic nucleon 
DA is unable to describe the nucleon form factors On the other hand, asymmetric DAs 
constructed via moments determined by local QCD sum rules following ^l|, as, for example, in 
4ll |42j (see on the left part of Fig. 0] for an illustration), seem to yield to strongly 



suppressed results for the magnetic nucleon form factor when transverse momentum — intrinsic 
and Sudakov — effects are included l36l |43| . Valuable information on the inner structure of the 
nucleon was recently provided in |44| in the context of instantons, where it was shown that 
the shape of the proton DA is far from the asymptotic one. 

While the nonperturbative nature of the nucleon is yet not well-understood, its evolution on 
the basis of the renormalization-group equation can be performed to a high level of accuracy 
within QCD perturbation theory. Indeed, within the basis of symmetrized Appell polynomials 



36l I45| , the nucleon evolution equation can be solved by employing factorization of the depen- 
dence on the longitudinal momentum from that on the external (large) momentum scale Q 2 
up to any desired polynomial order. 2 The spectrum of the corresponding anomalous dimen- 
sions of trilinear quark operators was also determined |36l . 45l 46l | and its large-order behavior 



seems to increase logarithmically, reflecting the enhanced emission of soft gluons that forces 
the probability for finding bare quarks to decrease (see Fig. HI right panel). This spectrum 
can be reproduced by the logarithmic fit 

7n (M) = c + d\n(M + b) . (4.1) 

The upper envelope of the spectrum is best described by the following values of the parameters 
with their errors: b = 1.90989±0.00676, c = -0.637947±0.000634, and d = 0.88822±0.000119. 
For the lower envelope, the corresponding values are b = 3.006 ± 0.483, c = —0.3954 ± 0.0290, 
and d = 0.59691 ± 0.00545. The spacing of eigenvalues at very large order is reproduced 
by the values b = -0.027 ± 0.728, c = -0.2460 ± 0.0248, and d = 0.291883 ± 0.00475. For 



2 The eigenfunctions of the nucleon evolution equation are linear combinations of symmetrized Appell poly- 
nomials, appropriately orthonormalized 
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FIG. 4: Left: The heterotic nucleon distribution amplitude, proposed in |41|. Right: Spectrum of the 
anomalous dimensions of trilinear twist-3 quark operators up to order M = 400. The solid lines (upper and 
lower envelopes of the spectrum) represent logarithmic fits up to the maximum considered order 400, taking 
into consideration all orders above 10. The dashed line gives for comparison a previous logarithmic fit |3Jj 
which takes into account all orders up to 150. 



every order M, there are M + 1 eigenfunctions of the same order with an excess of symmetric 
(under the permutation P 13 ) terms (denoted by black dots in Fig. 0]) by one for even orders. 
The total number of eigenfunctions up to order M is n max (M) = \[M + 1)(M + 2) and the 
corresponding (M + 1) eigenvalues are obtained by diagonalizing the (M + 1) x (M + 1) 
matrix. Up to order 150, both sectors (corresponding to the permutation parity S n = ±1) 
of eigenvalues are included. Beyond that order, for reasons of technical convenience, only the 
antisymmetric (open circles) ones have been taken into account. The multiplet structure of 
the anomalous dimensions spectrum was found independently later on in the context of a 
Hamiltonian approach to the one-dimensional XXX Heisenberg spin magnet of non-compact 
spin s = —1. 



V. CONCLUDING REMARKS 

Our discussion of the pion DA in the context of QCD sum rules with non-local condensates 
shows that the vacuum non-locality parameter can serve to extract valuable information on 
the underlying nonperturbative dynamics. The double-humped shape with suppressed end- 
points of the derived pion DA is in good agreement with the CLEO data with a la accuracy 
and agrees with the CELLO data as well. Progress of the non-local sum-rules approach to 
encompass tree-quarks states, like the nucleon, appears promising, while the perturbative 
apparatus for the evolution of such DAs is already well-developed. 

In conclusion, let us mention as a personal statement that the major part of our scientific 
work depends to a great extent on the power of factorization theorems and their usage in 
QCD in the context of form factors, structure functions, etc. Therefore, we feel particularly 
attached to Prof. Efremov, given also that he was the Leader of the BLTPh QCD group, where 
two of us (A.P.B. and S.V.M.) have been working for over a decade, and he was also one of the 
opponents of one of us (N.G.S.) in defending his Doctor fiziko-matematicheskih nauk degree. 
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